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We review the present status of CP violating problem in supersymmetric extensions of 
the standard model. We analyze the constraints imposed by the experimental limits of 
the electron, neutron, and mercury electric dipole moments on the supersymmetric CP 
phases and show that only the scenarios with flavour-off- diagonal CP violation remain 
attractive. These scenarios require hermitian Yukawa matrices which naturally arise in 
models with left— right symmetry or a 5(7(3) flavour symmetry. In this case, ek and e'/e 
can be saturated by a small non— universality of the soft scalar masses through the gluino 
and chargino contributions respectively. The model also predicts a strong correlation 
between Acp{b — > 57) and the neutron electric dipole moment. In this framework, the 
standard model gives a the leading contribution to the CP asymmetry in U — > "tjiKs 
decay, while the dominant chargino contribution to this asymmetry is < 0.2. Thus, no 
constraint is set on the non-universality of this model by the recent BaBar and Belle 
measurements. 



1. Introduction 

The understanding of the origin of CP violation is one of the remanining open ques- 
tions in particle physics. In the standard model (SM), CP violation and flavour tran- 
sition arise from the complex Yukawa couplings which have a Cabibbo-Kobayashi- 
Maskawa (CKM) mixing matrix with physical Sckm phase of order unity. At 
present, CP violation is observed in the K and B systems and the experimen- 
tal results are consistent with the SM. However, CP violation can be tested from 
the existing bounds on the electric dipole moment (EDM) of the neutron and the 
electron. It is remarkable that the SM contribution to the EDM of the neutron 
is of order 10~^*^ e.cm, while the experimental limit is of order 10^^^ e cm U. We 
expect that with the further improvements of experimental precision the EDM will 
provide a crucial test of CP violation. Also more experimental information on the 
CP violation will be obtained from the B-factory soon, which would be a crucial 
test of the CP violation in SM. Furthermore, SM model is unable to explain the 
cosmological baryon asymmetry of our universe, the presence of new sources of CP 
violation is required for this explanation. 

In supersymmetric (SUSY) extensions of the SM there are additional sources of 
CP violation, due to the presence of new CP violating phases which arise from the 
complexity of the soft SUSY breaking terms and the SUSY preserving /i-parameter. 
These new phases have significant implications and can modify the SM predictions 
in CP violating phenomena. In particular, they would give large contributions to 
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the electric dipole moment (EDM) of the electron, neutron and mercury atom 13, to 
CP violating parameters {ek, s'/^) of K ~ K system 0, and to the CP asymmetries 
in the B — B system □. These phases can be classified into two categories. The first 
category includes flavor-independent phases such as the phases of the /i-parameter, 
i3-parameter, gaugino masses and the overall phase of the trilinear couplings. The 
other category includes the flavor-dependent phases, i.e. the phases of the off- 
diagonal elements of Aij after the overall phase is factored out and phases in the 
squark mass matrix rafy Two of the flavor-independent phases can be eliminated 
by the U{1)r and U{1)pq transformations under which these parameters behave 
as spurious. The Peccei-Quinn transformation act on the Higgs doublets and the 
right-handed superfields in such a way that all the interactions but which mix the 
two doublets are invariant. The Peccei-Quinn charges are QpgifJ-) = Qpq{B^), 
Qa{i^) — Qpgi^rii) = 0. The U{l)p trnasforms the Grassmann variable — > 0e'" 
and the fields in such a way that the integral of the superpotential over the Grass- 
mann variables is invariant, i.e., the J7(l)jj charge of the superpotential is 2. As a 
result, Qr{B^) = Qpiy) ~ 2, Qr{A) = QR{mi) = —2. The six physical CP-phases 
of the theory are invariant under both U{1)r and U{1)pq, and can be chosen as 

Arg(A>,), Arg((i3M)*Ma), (1) 

where i = 1,2,3 and a ~ d,u,l. All other CP-phases can be expressed as linear 
combinations. If the A-terms and the gaugino masses are universal, there are two 
physical phases Arg(A*m), Arg{B*A). 

However, the non-observation of EDMs imposes a stringent constraint on flavor- 
independent SUSY phases. Putting these new phases to zero is not natural in the 
sense of 't Hooft u since the Lagrangian does not acquire any new symmetry in the 
limit where these new phases vanish. This is the so-called SUSY CP-problem. In this 
article we will review the constraints imposed by the EDMs on the flavor-diagonal 
CP-phases and the possible scenarios allowing to supress the EDM contributions: 
small SUSY CP phases, heavy sfermions, EDM cancellation, and flavor off-diagonal 
CP violation. Also we will discuss aspects of CP violation in the K and B systems 
due to the flavor-off-diagonal phases from non-degenerate A-terms. 

This article is organized in the following way. In section 2 we briefly discuss the 
origin of the soft SUSY breaking terms. Section 3 is devoted to the study of the 
constraints from the EDMs of the electron, neutron, and mercury atom in generic 
SUSY models. We present the possible ways to avoid overproduction of EDMs and 
explain that, in SUSY models with flavor off-diagonal phases, the EDMs can be kept 
sufficiently small while these phases unconstrained. The effect of the off-diagonal 
phases in the CP violation process in the kaon system is given in section 4. In section 
5 we analyze the CP violation in the i3-sector and show that in this framework 
the large CP asymmetry in the B ipKs decay is given by the SM contribution 
while the SUSY contribution is very small. In contrast, the SUSY contribution to 
the CP asymmetry in the B Xsj decay can be as large as ±10%. Finally, the 
conclusions are presented in section 6. 
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2. The origin of the soft terms 

It is clear that the SUSY CP violating problem is a problem of SUSY breaking since 
the relevant phases originate from SUSY breaking terms. In this section we briefly 
discuss the possible mechanisms which may give rise to the SUSY soft breaking 
terms. The general structure one has in mind includes three sectors: 

i) The observable sector which comprises all the ordinary particles and their SUSY 
partners, 

ii) a "hidden" or "secluded" sector where the breaking of SUSY occurs, 

iii) the messengers of the SUSY breaking from the hidden to observable sector. 
The two most explored alternatives that have been studied in this context are: 

a) SUSY breaking supergravity where the mediators are gravitational interactions, 
the scale of SUSY breaking is Ms — \/Mp Mw and the mass of the fermionic 
partner of the gravition, the gravitino, is m^^2 Mg/Mp. 

b) SUSY broken at a much lower scale with messengers provided by some gauge 
interactions. In this case 7713/2 = Mg/Mp is generally very small, the gravitino 
being the most likely candidate for the LSP of the theory. 

We start the former and still more popular alternative. As well known that 
the invariance under local SUSY transformation implies invariance under local co- 
ordinate change. Thus, local SUSY (supergravity) naturally includes gravity. The 
effective action, up to two derivatives, is completely specified in terms of three 
functions which depend on the chiral superfields of the theory. The Kahler poten- 
tial K{^,^^), the analytic superpotential H^(<I>) and the gauge kinetic functions 
/"(<!>) which are also analytic functions of the chiral superfields. The tree level 
supergravity scalar potential is given by □ 



G^{G-')lG,-3\+\fjD"D0, (2) 

where G($,$*) = + log|VF($)|2, G' = and = g^§^. The grav- 

itino mass is given by Eqjl shows that unlike the case of the global supersymmetry 
one can simultaneously have (Gi) ^ ^ (which breaks SUSY) and vanishing cos- 
mological constant. 

Supergravity is broken in a "hidden sector", namely a sector of the theory 
that couples to the "observable sector" of quarks, leptons, gauge fields, Higgs 
and their supersymmetric partners, only through gravitational interactions. The 
simplest model of a supersymmetry breaking hidden sector is known as "Polonyi 
model". In this model the minimal form of the Kahler potential K = z*^Zi + y*'^yr is 
assumed {zi and denote the fields of hidden and observable sectors respectively). 
The resulting soft breaking terms are obtained taking the so-called flat limit , i.e. 
sending Mp 00 while keeping the ratio Mg/Mp = TO3/2 fixed ( Ms denotes the 
scale of supergravity breaking) . One obtains a common mass term for all the scalar 



^ The auxiliary fields of the chiral and vector supermultiplets are given in terms of d see Ref.t 
One can see that if at least one of the d VEV's is non-vanishing SUSY is broken 
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particles in the observable sector which is equal to the gravitino mass. Also in this 
model we get a universal A term {i.e. Ajj = Ad = Al = A). 

However, it is possible to obtain effective potentials in which this universality 
is absent by taking the kinetic terms for the chiral superfields to non-minimal. As 
recently stresse, the soft supesymmetry breaking parameters may be non-universal 
in the effective theories which are derived from the superstring theories. In general 
supergravity, the soft scalar mass and the v4-parameter are given by 

K = ^1/2 - F"'F^d^dn In , (3) 

where m refers to the SUSY breaking fields and = K^a- Note that in Eq.(^ 
we have diagonal soft scalar masses due to the assumption that the Kahler metric 
is diagonal. General supergravity models allow for a non-vanishing off-diagonal 
Kahler metric Kap- However, such mixing typically does not appear in superstring 
models at the leading order due to additional "stringy" symmetries beyond those 
of the Standard Model. The gaugino masses are given in terms of the gauge kinetic 
function fa, where Re/a = and the subscript a represents the corresponding 
gauge group. 

Ma = ^(Re/a)-'i^"5„/, . (5) 

Apart from ma, Aafj^, and Ma soft parameters, soft bosonic bilinear can be also 
present. In case of the MSSM, all the symmetries of the low energy allow for a 
superpotential term couphng the two Higgs doublets of the form W = fiHiH2. The 
associated soft breaking term in the scalar potential will have the form B11H1H2 + 
h.c, where _B is a dimensionful parameter of order gravitino mass. It is well known 
that in order to get appropriate SU{2)l x U{1)Y breaking, the has to be of the 
same order of magnitude as the SUSY breaking soft terms. 

Now, we turn to the alternative mechanism of breaking supersymmetry at low 
energies and that gauge interactions are the "messenger" of supersymmetry break- 
ing lI. The gauge mediated supersymmetry breaking (GMSB) models have several 
attractive features. As gauge interactions are flavor blind, squark and slepton masses 
are universal. In the minimal model of this kind, the messanger fields transform as 
a single flavor of 5 -I- 5 of SU{5). Hence, there are SU{2)l doublets I and I and 
SU (3)c triplets q and q. In order to introduce supersymetry breaking into the mes- 
sanger sector, these flelds may be coupled to a gauge singlet spurion, S, through 
the superpotential 

W = X1SII + \2Sqq, (6) 



due to interactions with some supersymmetry breaking sector of theory. The field 
S has a non-zero expectation value both for its scalar and auxiliary component, {S) 
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and (Fs). Integrating out the messenger sector give rise to gaugino masses at one 
loop. 

"^A. -Q^A, (7) 

where A — ^j^, ci ~ 5/3, C2 = C3 = 1 and ai — ^^Jig^^, ■ For the scalar masses one 
has 



(8) 



where C3 = 4/3 for color triplets and zero for singlets, C2 = 3/4 for weak doublets 
and zero for singlets, and Y — 2(Q — Ts) is the ordinary hypercharge. In this model, 
the v4-terms are high order effects. The /i and B parameters can arise through 
interactions of the Higgs field with various singlets El. This model has few new 
parameters relative to the minimal standard model. These parameters are A, /i and 
B. Since the scalar masses are functions of only gauge quantum numbers, these 
models also automatically solve the supersymmetric flavor problem. 

The most important difference between the models of the gravity mediated 
SUSY breaking and the models of gauge mediation is that in the later models the 
LSP is the gravitino and its mass is of order 10~^° GeV. Thus leads to a very 
different phenomenology from the formare model. It allows the lightest neutralino 
Xi to decay to photon plus gravitino and this gives a signal for SUSY as a missing 
energy. 



3. SUSY CP violating phases and constraints from the EDMs 

As mentioned in the introduction, in supersymmetric theories there are several 
possible sources for CP violating phases in addition to the CKM phase. These CP 
violating phases have an important impact on the phenomenology of CP violatin. 
They can induce an EDM of quarks andjeptons at one-loop level. The EDMs of 
the neutron, electron, and mercury atom El impose severe constraints on the flavor- 
diagonal phases, the so-called SUSY CP problem. 

dn < 6.3 X 10"^^ e cm , (9) 
de < 4.3 X 10^2^ e cm , (10) 
dHg < 2.1 X 10^2* e cm . (11) 

With the expected improvements in experimental precision, the EDM is likely to 
be one of the most important tests for physics beyond the Standard Model for some 
time to come, and EDMs will remain a difficult hurdle for supersymmetric theories. 
Indeed it is remarkable that the SM contribution to the EDM of the neutron is of 
order 10~^° e cm, whereas the "generic" supersymmetric value is 10~^^e cm. 

In this section we review EDM constraints in the context of supersymmetric 
theories as well as known mechanisms to suppress the EDMs. We also study to 
what extent different scenarios rely on assumptions about the neutron structure, 
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i.e. chiral quark model vs parton model. Let us first summarize the contributions 
to the three most significant EDMs, beginning with the most reliable, the electron 
EDM. 

The electron EDM is defined by the effective CP-violating interaction 

C = -^4ea,,,75e F^'' , (12) 

where -F'"^ is the electromagnetic field strength. The experimental bound on the 
electron EDM is derived from the electric dipole moment of the thallium atom and 
is given by 

4 < 4 X 10"27(, ^ ^i^-j 

In supersymmetric models, the electron EDM arises due to CP-violating 1-loop 
diagrams with the chargino and neutralino exchange: 

4 = df + . (14) 

Since the EEDM calculation involves little uncertainty it allows to extract reliable 
bounds on the CP-violating SUSY phases. 

The neutron EDM has contributions from a number of CP-violating operators 
involving quarks, gluons, and photons. The most important ones include the electric 
and chromoelectric dipole operators, and the Weinberg three-gluon operator: 



£ = - -d^qa^.J^q F^" - -d'^^ qa^.-f^T^q G^^^ 



'^d^qcT^^.lsq F^" - '-i 

" -^dP ] abcG aiipG^^jG cXa^^^^" , (15) 

where Ga^i/ is the gluon field strength, and fabc are the SU(3) generators and 
group structure coefficients, respectively. Given these operators, it is however a 
nontrivial task to evaluate the neutron EDM since assumptions about the neutron 
internal structure are necessary. In what follows we will study two models, namely 
the quark chiral model and the quark parton model. Neither of these models is suffi- 
ciently reliable by itself u, however a power of the combined analysis should provide 
an insight into implications of the bound on the neutron EDM. AnotJier approach 
to theneiitron EDM based on the QCD sum rules has appeared in lid and earlier 
work llil, We note that in any case the NEDM calculations involve uncertain 
hadronic parameters such as the quark masses and thus these calculations have a 
status of estimates. The major conclusions of the present work are independent of 
the specifics of the neutron model. 

The chiral quark model is a nonrelativistic model which relates the neutron 
EDM to the EDMs of the valence quarks with the help of the SU(6) coefficients: 

4 1 

dn = -^dd - -du . (16) 
The quark EDMs can be estimated via Naive Dimensional Analysis Ei as 



d,^,-d- + ,-f^d- + ,-'-^d-, (17) 



A 
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where the QCD correction factors are given by rj ~ 1.53, r;*-^ ~ 77"-' ~ 3.4, and 
A ~ 1.19 GeV is the chiral symmetry breaking scale. We use the numerical val- 
ues for these coefficients as given in EJ. The parameters rp^'^ involve considerable 
uncertainties steming from the fact that the strong coupling constant at low ener- 
gies is unknown. Another weak side of the model is that it neglects the sea quark 
contributions which play an important role in the nucleon spin structure. 

The supersymmetric contributions to the dipole moments of the individual 
quarks result from the 1-loop gluino, chargino, neutralino exchange diagrams 

^E.C ^ iE,C) ^ (£,C) ^ iE,C) ^ (^g) 

and from the 2-loop gluino-quark-squark diagrams which generate . 

The parton quark model is based on the isospin ^mmetry and known contri- 
butions of different quarks to the spin of the proton cB. The quantities A, defined 
as {n\^qjf_i'j5q\n) = Aq S*^, where S'^ is the neutron spin, are related by the isospin 
symmetry to the quantities (Ag)p which are measured in the deep inelastic scat- 
tering (and other) experiments, i.e. A„ = (Ajj)p, A^j = (Ai,)p, and A^ = (As)p. To 
be exact, the neutron EDM depends on the (yet unknown) tensor charges rather 
than these axial charges. The main assumption of the model is that the quark 
contributions to the NEDM are weighted by the same factors Ai, i.e. Ej 

d„ =7y^(Addf + A„df -f Asdf) . (19) 

In our numerical analysis we use the following values for these quantities = 
0.746, Au = -0.508, and A^ = -0.226 as they appear in the analysis of RefJla. As 
before, we have 

= 4 (^) + df + df . (20) 

The major difference from the chiral quark model is a large strange quark contribu- 
tion (which is likely to be an overestimate i) . In particular, due to the large strange 
and charm quark masses, the strange quark contribution dominates in most regions 
of the parameter space. This leads to considerable numerical differences between 
the predictions of the two models. 

The EDM of the mercury atom results mostly from T-odd nuclear forces in 
the mercury nucleus, which induce the effective interaction of the type (I • V)(5(r) 
between the electron and the nucleus of spin I In turn, the T-odd nuclear Jorces 
arise due to the effective 4-fermion interaction ppnij^n. It has been argued S that 
the mercury EDM is primarily sensitive to the chromoelectric dipole moments of 
the quarks and the limit 

dHg < 2.1 X IQ-^^e cm (21) 

can be translated into 

\dd - du ~ 0.012<if |/5s < 7 X IQ-^^cm , (22) 

where gs is the strong coupling constant. As in the parton neutron model, there 
is a considerable strange quark contribution. The relative coefficients of the quark 
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contributions in (|2^) arc known better than those for the neutron, however the 
overaU normahzation is still not free of uncertainties ll. 

Below we list formulae for individual supersymmetric contributions to the EDMs 
due to the FeynmaiL diagrams in Fig.^ In our presentation we follow the work of 
Ibrahim and Nath llil. 




Fig. 1. Leading SUSY contributions to the EDMs. The photon and gluon lines are to be attached 
to the loop in all possible ways. 



Neglecting the flavor mixing, the electromagnetic contributions to the ferniion 
EDMs are given by0: 



7X+ (E) 



/e 



fc=i 



An sin^ 9w 



47r sin^ 6w 



5^5^Im(r™fe) 



fc=l i=l 
2 2 



5^5^Im(r,,fe) 



fc=l i=l 



Ml 



Wk 



Qu B 



dk 



Ml 



{Qu - Qa) A 
[Qd - Qu) A 



dk 



Ml. 



An sin^ 6 



V ^Im(re,) A 



Here 



.J2J2MVf^k)^QfB 



fk 



(23) 



q2k^qlk I 



(24) 



with 03 being the gluino phase and Dg defined by DlM?Dq = diag(M|^ , M|J. The 



sfermion mass matrix M? is given by 
M? = (Ml^+'^ f^+MUh^Qf sir? Ow) cos 2/3 



M'l.+mf'^ + M'^QfSm^ew cos2/3 



,(25) 
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where Rf — cot /3 (tan/?) for I3 — 1/2 (—1/2). The chargino vertex T fik is defined 
as 



(26) 



and analogously for the electron; here U and V are the unitary matrices diagonal- 
izing the chargino mass matrix: U* M^+V~^ = diag(m^+, m^+). The quantities k/ 
are the Yukawa couplings 



V^iTT-w sin (3 ' 
The neutralino vertex 77/,;^ is given by 



md,e 



V2mw cos (3 



(27) 



Vfik 



~V2{taiiew{Qf - h,)Xu + hjX2^}D}^k - ^^fXuD}2k 
V2ta,ii6wQfXuDf2k - HfX^Dfik 



(28) 



where is the third component of the isospin, 5 = 3 (4) for — —1/2 (1/2), and 
X is the unitary matrix diagonalizing the neutralino mass matrix: M^oX — 
diag(m^o , m^o , m^o , m^o ) . In our convention the mass matrix eigenvalues are posi- 
tive and ordered as m^o > to^o > ... (this holds for all mass matrices in the paper). 
The loop functions A{r),B{r), and C{r) are defined by 

21nr' 



A{r) = 
B(r) = 
C{r) = 



1 

2(1 -r)2 
1 

2(r-l)2 
1 

6(r- 1)2 



3-r 
1 + r 



1 - r 
2r Inr 



lOr - 26 



1 - r 
2r In r 



1 



r 



18 Inr 

1 -r 



(29) 



The chromoelectric contributions to the quark EDMs are given by 

iC) _ 9sas 



-9 9. 
167r2 



k=l 
2 2 2 



9fc 



c 



Mi 



M? 



EEMr,.fe)-r^B 



fe=i j=i 



qk 



M? 
qk 



fe=i j=i 



Af,2 

qk 



B 



to2o 
Xi_ 

M? 

qk 



(30) 



Finally, the contribution to the Weinberg operator from the two-loop gluino- 
top-stop and gluino-bottom-sbottom diagrams reads 

3 



' 47r 



Im(ri2) ^L_^ H(zi,Z2,zt) + (t^b), (31) 
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where Zi 



M3 



M3 



The two-loop function H(zi, Z2, zt) is given by 



H{zi,Z2,zt) = - f dx f du f dy x{l 
^ Jo Jo Jo 



■ x)u 



NIN2 

£)4 



(32) 



where 



Ni = u(l — x) + Ztx{l — x){l — u) — 2ux[ziy + 22(1 — y)] 



No 



{l-xf{l-tif + u''-^x\l-uf 
9 



(33) 



D — u{\ — x) + Ztx{l — x){\ — u) + ux[ziy + 2:2(1 — y)] . 

The numerical behaviour of this function was studied in 111. We emphasize that the 
b-quark contribution is significant and often exceeds the top one. 

In addition to the Weinberg two-loop diagram, there is another (liajr-Zee) two- 
loop contribution which originaies from the CP-odd Higgs exchange c2l. Its numer- 
ical effect is however negligible □. 

Let us now discuss possible ways to avoid overproduction of EDMs in super- 
symmetric and string-inspired models. There are four known scenarios allowing to 
suppress the EDM contributions: small SUSY CP-phases, heavy sfermions, EDM 
cancellations, and flavor-off-diagonal CP violation. 



1 Small SUSY CP-phases: 

For a light (below 1 TeV) supersymmetric spectrum, the flavor-independent 
SUSY CP phases have to be small in order to satisfy the experimental EDM 
bounds □. 





Fig. 2. EDMs as a function of tpA (left) and (right). 1 - electron, 2 - neutron (chiral model), 
3 - mercury, 4 - neutron (parton model). The experimental limit is given by the horizontal line. 
Here tan/3 = 3, mo = mi/2 = A = 200 GeV. 
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In Figsj^ we illustrate the EDMs behaviour as a function of the CP-phases in 
the mSUGRA-type models, where we have set mo = mi/2 = A = 200 GeV. At low 
tan/3, the EDM constraints impose the following bounds (at the GUT scale): 



< 10-2 - IQ-i , 

(I^AU < 10"' • (34) 

We note that 4>a is less constrained than c/)^ and (pMt- There are two reasons for 
that: first, (j)A is reduced by the RG running from the GUT scale down to the 
electroweak scale and, second, the phase of the {Sfj^)LR mass insertion which gives 
the dominant contribution to the EDMs is more sensitive to (j)^ and (j)Mi due to 
1^1 < /itan/3. We note that the bounds ( ^4|) stay practically the same if we allow 
for an "uncertainty" factor 2-3 in the overall EDM normalization. 

Generally it is quite difficult to explain why the soft CP-phases have to be small. 
In principk, small CP-phases could appear if CP were an approximate symmetry 
of nature E2l. However ,_Lecent experimental results show that CP violation in the 
B — B mixing is large I2I and thus the approximate CP hypothesis cannot be mo- 
tivated. 



2- Heavy sfermions: 



This possibility is based on the decoupling of heavy supersymmetric particles. 



Even if one allows 0(1) CP violating 
SUSY spectrum is sufficiently heavy I 



rases, their effect will be negligible if the 
Generally, SUSY fermions are required 



to be lighter than the SUSY scalars by, for example, cosmological considerations. 
So the decoupling scenario can be implemented with heavy sfermions only. Here 
the SUSY contributions to the EDMs are suppressed even with maximal SUSY 
phases because the sauarks in the loop are very heavy and the mixing angles are 
small. It was shown □ that in order to satisfy the EDM of the mercury atom the 
sfermion masses have to be of order 10 TeV, which leads to a large hierarchy between 
SUSY and electroweak scales. In Fig.^ we display the EDMs as functions of the 
universal scalar mass parameter mg for the mSUGRA model with maximal CP- 
phases 0^ = = 7i'/2 and mi/2 = A = 200 GeV. We observe that all EDM 
constraints except for that of the electron require mo to be around 5 TeV or more. 
The mercury constraint is the strongest one and requires mo ~ 10 TeV. 

This leads to a serious fine-tuning problem. Recall that one of the primary mo- 
tivations for supersymmetry was a solution to the naturalness problem. Certainly 
this motivation will be entirely lost if a SUSY model reintroduces the same problem 
in a different sector, i.e., for example a large hierarchy between the scalar mass and 
the electroweak scale. Also note that the dependence of the decoupling scale on the 
EDMs is quite slow. For instance, if we relax the mercury bound by a factor of 2 
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Fig. 3. EDMs as a function of the universal mass parameter rriQ. 1 - electron, 2 — neutron (chiral 
model) , 3 - neutron (parton model) , 4 — mercury. The experimental limit is given by the horizontal 
line. Here tan^ = 3, mi/2 = A = 200 GeV, 0^ = (/)^ = 7r/2. 

the decoupling scale changes from 10 TeV to about 8 TeV. Thus our conclusions 
are insensitive to the exact values of the EDMs. 

3- EDM cancellations: 

The cancellation scenario is based on the fact that large cancellations among 
different^amtiibutions to the EDMs are possible in certain regions of the parameter 
space oHlHil'E^ which allow for 0(1) flavour — independent CP phases. However, a 
recent analysis 13 has shown that this possibility is practically ruled out if in addition 
to the electron and neutron EDM constraints, the mercury constraint is imposed. 

For the case of the electron, the EDM cancellations occur between the chargino 
and the neutralino contributions. For the case of the neutron and mercury, there 
are cancellations between the gluino and the chargino contributions as well as can- 
cellations among contributions of different quarks to the total EDM. It has been 
shown i that the parameters allowing the EDM cancellations strongly depend on 
the neutron model. For example, in the parton model, it is more difficult to achieve 
these cancellations due to the large strange quark contribution. Therefore, one can- 
not restrict the parameter space in a model-independent way and caution is needed 
when dealing with the parameters allowed by the cancellations. 

In mSUGRA, the EDM cancellations can occur simultaneously for the electron, 
neutron, and mercury along a band in the {4>A, (t>n) plane (Fig.^. However, in this 
case the mercury constraint requires the /i nhase to be 0(10^^) and the magnitude 
of the A-terms to be suppressed (~ O.lmo) a which results in only a small effect of 
the A-terms on the phase of the corresponding mass insertion. Thus the fact that the 
phase of the A-terms is unrestricted should not be attributed to the cancellations, 
but rather to the suppressed contribution of the A-terms. 

If the gluino phase is turned on, simultaneous EEDM, NEDM, and mercury 
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Fig. 4. Phases allowed by simultaneous electron, neutron, and mercury EDM cancellations in 
mSUGRA. The chiral quark neutron model is assumed. Here tan/3 = 3, mo = mi/2 = 200 GeV, 
A = 40 GeV. 



7t/10 

% 

Fig. 5. Bands allowed by the electron (red), neutron (green), and mercury (blue) EDMs can- 
cellations in the mSUGRA-type model with nonzero gluino and bino phases. Here tan (3 = 3, 
mo = mi/2 = 200 GeV, A = 40 GeV, (f>i = (jj^ = tt/10. 

EDM cancellations are not possible. The gluino phase affects the NEDM cancella- 
tion band while leaving the EEDM cancellation band almost intact. An introduction 
of the bino phase (f>i qualitatively has the same "off-setting" effect on the EEDM 
cancellation band as the gluino phase does on that of the NEDM .(Fig.||). Note 
that the bino phase has no significant effect on the neutron and mercury cancel- 
lation bands since the neutralino contribution in both cases is small. When both 
the gluino and bino phases are present (and fixed), simultaneous electron, neutron, 
and mercury EDMs cancellations do not appear to be possible along a band. This 
conclusion has been also obtained in type I string inspired models 

The cancellation scenario involves a significant fine-tuning. Indeed, restricting 
the phases to the band where the cancellations occur does not increase the symme- 
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try of the model and thus is unnatural according to the t'Hooft's criterion i. 
4 Flavor- off- diagonal CP violation: 

Finally, we consider the possiblity that the SUSY CP violation have a flavor off- 
diagonal character just as in the SM cB. In the SM, CP violation appears in flavor- 
changing processes which is one of the reasons why the predicted EDMs are small. A 
similar situation may occur in supersymmetric models. In this case, the origin of the 
CP violation is assumed to be closely related to the origin of flavor structures rather 
than the origin of SUSY breaking. Thus the flavor blind quantities as the /i-term 
and gaugino masses are real. As mentioned, the EDMs of the neutron, electron, 
and mercury atom impose severe constraints on the flavor-diagonal phases. Indeed, 
these are the only phases that enter into the EDM calculations. However, this leaves 
the possibility that the flavor-off-diagonal phases are non-zero. Of course, one has 
to specify the basis in which the considered quantity is classified as off-diagonal. 

This class of models requires hermitian Yukawa matrices and A-terms Ej, which 
forces the flavor-diagonal phases to vanish (up to small RG corrections) in any 
basis. To see this, first note that the trilinear coupings 

A?, = A%Y[; (35) 

are also hermitian. The quark Yukawa matrices are diagonalized by the unitary 
transformation 

q^Vq , Y"^ {V^f {V^y = diag(/i«,/i«,/i|) , (36) 

such that the CKM matrix is Vckm — {V^yV^. If we transform the squark fields 
in the same manner, which defines the super-CKM basis, we will have 

i« (Y'lf i« (F«)* . (37) 

As a result, the trilinear couplings remain hermitian and the flavor-diagonal CP- 
phases inducing the EDMs vanish. Note that this argument would not work if in 
the original basis the trilinear couplings but not the Yukawas were hermitian since 
the diagonalization would require a biunitary transformation which would generally 
introduce the diagonal phases. 

The hermiticity is generally spoiled by the renormalization group (RC) running 
from the high energy scale down to the electroweak scale (see, e.g. E3). In models 
under consideration we have the following setting at high energies: 

Arg(Mfc) = Arg(B) = Argifi) = . (38) 

Generally, the off-diagonal elements of the A-terms can have 0(1) phases without 
violating the EDM constraints. Due to the RG effects, large phases in the soft tri- 
linear couplings involving the third generation generate small phases in the flavour- 
diagonal mass insertions for the light generations, and thus induce the EDMs. How- 
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ever, the amount of generated "non-hermiticity" depends on A'^,A'^ 
etc. and thus is suppressed by the ofF-diagonal elements of the Yukawa matrices. 
In particular, the size of the flavor-diagonal phases is quite sensitive to Y13. For 
^13 ^ 0(10-3), the RG effects are unimportant and the induced EDMs are below 
the experimental limits. 

In supergravity models, the trilinear parameters are given in terms of the Kahler 
potential and the Yukawa couplings 

F"" [k^ + In Y^i3^ - Hk^kpk^)] , (39) 

where the Latin indices refer to the hidden sector fields that break SUSY and the 
Greek indices refer to the observable fields. According to our previous assumption, 
the is real. Also k^ and kj^ are always real, thus the A-terms are Hermitian 
if the derivatives of the Kahler potential are either generation-independent or the 
same for the left and right fields of the same generation, i.e., if Kqj^ , Kjjj^ , — 

kqz, Kur. ■ These conditions are usually satisfied in string models. It is interesting 
to note that although the SUSY breaking does not bring in new source of CP 
violating, the trilinear soft parameters involve off-diagonal CP violating phases of 
0(1). This stems from the contribution of the term dm^nYap^, whiciLhas been 
found to be significant and sometimes even dominant in string models E3. In what 
follows, we will show that these phases are unconstrained by the EDMs and will 
study their phenomenological implications in the K and B systems. 

The relevant quantities appearing in the soft Lagrangian are iYi^)ij = {Yq)ij {Aq)ij 
(indices not summed) which are also Hermitian at the GUT scale. For the sake of 
definiteness, we consider the following Hermitian Yukawa matrices at the GUT scale 

/ 5.94x10-'' 10-3 i -2.03 X 10-2 \ 
r" = A„ -10-3 i 5.07 X 10-3 2.03 x IQ-^ i , (40) 
V-2.03 X 10-2 -2.03 X 10-5 i 1 / 

/ 6.84 x 10-3 (1.05 + 0.947 i) X 10-2 -0.023 \ 
= Ad (1.05 - 0.947 i) x IO-2 0.0489 0.0368 i , (41) 

V -0.023 -0.0368 i 1 / 

where A„ — mt /v sin P and A^ = mt,/u cos /3. These matrices reproduce, at low 
energy, the quark masses and the CKM matrix. The renormalization group (RG) 
evolution of Yukawa couplings and the A terms slightly violate the Hermiticity. 
Therefore, the resulting Y,^ at the electroweak scale has very small non-zero phases 
in the diagonal elements (due to the large suppression from the off-diagonal entries 
of the Yukawa). However, what matters is the relevant phases appearing in the 
squark mass insertions in the super-CKM basis, i.e., the basis where the Yukawa 
matrices are diagonalized by a unitary transformation of the quark superfields Ul,r 
and Dl h (Note that since the Yukawas are Hermitian matrices they are diagonal- 
ized by one unitary transformation, i.e. = ^r)- Accordingly the trilinear terms 
Y^ transform as Y^ V Y^V^ . Thus the Y,^ stay Hermitian to a very good 
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degree in the super-CKM basis. Therefore, the imaginary parts of the flavor con- 
serving mass insertions 



1 



d(u) 



that appear in the EDM calculations are suppressed. In the above formula the rriq 
refers to the average squark mass and Vi — (H^) /V^. The current experimental 
bound I leads to the constraint Im((5f j;"'')L_R 10^^. 

We start our analysis by revisiting the EDM constraints on the flavor off diagonal 
phases of SUSY models with Hermitian Yukawa as in Eq.(|4^) and the following 
Hermitian A-termaH3: 

/ An Ai2 e*"^!^ A13 e^'f^^ \ 
Ad = Au=[ A12 e-'^'^ A22 A23 e'^P^' . (43) 

We also assume that the soft scalar masses and gaugino masses Ma are given by 

Ma^mi/2, a = 1,2,3, (44) 

2222 /Ac:\ 

"^Q = rriH, = = m^, (45) 
mlj = mjj = ml diag{l, (5i, (52}, (46) 

where the parameters Si and Aij can vary in the ranges [0, 1] and [—3, 3] respectively. 
Note that in most string inspired models, the squark mass matrices are diagonal 
but not necessary universal. The non-universality of the squark masses is not con- 
strained by the EDMs. However, this non-universality (specially between the first 
two generations of the squark doublets) is severely constrained by AMk and Sk- 

In Fig. H we display scatter plots for the neutron and mercury EDMs versus 
the phase (pi2 for tan/3 = 5, mo = mi/2 = 200 GeV, Aij are scanned in the range 
[—3,3], and the phases ipis and (^23 are randomly selected in the range [0,7r]. As 
stated above, the parameters 5i are irrelevant for the EDM calculations and we 
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set them here to one. FinaUy, since /i is real the EDM results display very little 
dependence on tan/3. 

It is important to mention that we have also imposed the constraints which 
come from the requirement of absence of charge and colour breaking rninima as 
well as the requirement that the scalar potential be bounded from below cB. These 
conditions may be automatically satisfied in minimal SUSY models, however in 
models with non-universal yl-terms they have to be explicitly checked. In fact, 
sometimes these constraints are even stronger than the usual bounds set by the 
flavor changing neutral currents 

As can be seen from Fig. ^, the EDMs do not exceed the experimental bounds 
for most of the parameter space. Generally, they are one or two order of magni- 
tude below the present limit, and the flavor-off diagonal phases of the ^-terms can 
be 0(1) without fine-tunning. The points that lead to mercury EDM above the 
experimental bound correspond to ip23 — Tr/2. This phase induces a considerable 
contribution to the chromoelectric EDM of the strange quark C| . Thus the com- 
patibility with mercury EDM experiment requires that the phase (^23 should be 
slightly smaller than tt/2. 



4. CP violation in the Kaon system 

We have shown in the previous section that the Hermiticity of the Yukawa couplings 
and A-terms allows the existence of large off-diagonal SUSY CP violating phases 
while keeping the EDMs sufficiently small. However, the important question to 
address is whether these "EDM-free" phases can have any implication on other CP 
violation experiments. In this section, we will concentrate on possible effects in the 
kaon system. 



4.1. Indirect CP violation 

In the kaon system and due to a CP violation in K'^ — K'^ mixing, the neutral kaon 
mass eigenstates are superpositions of CP-even (Ks) and CP-odd (Kl) compo- 
nents. However, the CP-odd Kl decays into two pions through its small CP-even 
component. This decay, tttt, was the the first observation of CP violation. 

The measure for the indirect CP violation is defined as 

A{Kl ^ , . 

£k = -7777 r • (47) 

A[Ks ~* tttt) 

The experimental value for this parameter is Ek — 2.28 x 10~^. Generally, ek can 
be calculated via 



(48) 
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Here H^^='^ is the effective Hamiltonian for the AS' = 2 transition. It can be 
expressed via the operator product expansion as 

H^J^=^ = J2 C^iM + '^■c- ■ (49) 

i 

In the above formula, Ci(/i) are the Wilson coefficients and Qi are the EH local 
operators. 

The relevant operators for the gluino contribution are til 

Q4 = d%sldis% Qs = 4445fl, (50) 

as well as the operators Qi^2.3, that are obtained from (3i,2.3 by the exchange 
L i?. In the latter equations, a and /3 are SU{3)c colour indices, and the colour 
matrices obey the normalization Tr (f^t^) = 5°"^ /2. Due to the gaugino dominance 
in the chargino-smiark loop, the most significant contribution is associated with 
the operator Qi l3. 

In the presence of supersymmetric [g and x^) contributions, the following result 
for the amplitude M12 is obtained: 

Afi2 = Mf^ + + Mf^ . (51) 

The SM contribution can be written as US 

= ^{^f],mKBKT* , (52) 

where Bk is the renormalization group invariant which is given by 0.85 ± 0.15, 
Jk = 160 MeV is the K-meson decay constant, and rriK = 490 MeV the K-meson 
mass. The function !F can be give as 

T = TjiXlSaixc) + ri2XtSa{xt) + 2r]3XcXtSoixc, Xt) , (53) 

and = VH^d, Vz are QCD correction factors {rji = 1.38 ± 0.20, 7/2 = 0.57 ± 0.01, 
773 — 0.47 ± 0.04 - NLO values), and 5o(xi) are the loop functions given by 

(fji^ \ 1.52 
170GeV / 



.170GeV/ 

Xt 3xt _ 3xt log Xt 

4(1 -Xt) " Ail -xtY 



Soix,, Xt) = X, ( log ^ - -jji^ - ^ * (54) 



The supersymmetric term Aff2 is given by 
-as 1 



Mf2 - ^l£^3"^^/^i + (^12)Ir) Bi(p) (2AxMx)+66U{x) 



+ {5f,)LL{6i2)B.R 

+ B,i^i) (6xfe{x)+30Uix) 



Mk ^ ' 



B 



4(/i) (378a:/6(x)- 54/6(2:) 



2 , /.d N2 N ^ Mk 
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2 2 J \ms(n) + md(n) , 



-99B4(M)-45S5(/i) 



Mx)}, (55) 



where x = (rrig/mq) and the functions fe{x)^ feix) are given by 

- ^(^31)^ . (56) 

6n+3x)lnx + x^ -9x'^ -9x + 17 

- ^ 6(.-l)3 ■ 

The matrix elements of the operators Qi between the _fC-meson states in_the vacuum 
insertion approximation (VIA), where B ~ 1, can be found in Ref.E3. The VIA 
generally gives only a rough estimate, so other methods, e.g. lattice QCD, are 
required to obtain a more realistic value. The matrix elements of the renormalized 
operators can be written as E£l 

{K°\Q,i^i)\K°) = ^MKfKBiifi), (58) 

iK^m,)\K^) - ( ^.(,)?^,(,) )^^^^/i^3(.), (60) 

(if°|Q4(M)l^°) = ( , , MKfKB,{^,) , (61) 

where Qi{^i) are the operators renormalized at the scale /j. The expressions of the 
matrix elements of the operators Q1-3 are valid for the operators Qi-s and for 
fi^2 GeV we have @ 

Si(^)=0.60, S2(Ai)=0.66, 53(/i) = 1.05, 

S4(a*) = 1-03, BsIa^) 0.73 . 

Using these values, the gluino contribution to the K'^ — K'^ can be calculated via 
Eq.(|4^). As mentioned above, for universal soft scalar masses the LL and RR in- 
sertions are generated only through the RG running and can be neglected. The LR 
and RL mass insertions appear at the tree level and may have tangible effects. It is 
worth mentioning that, the RL and LR mass insertions contribute with the same 
sign in Eq.(|55|) and for Hermitian A-terms they are almost equal, so no cancellation 
between these two contributions occurs. 

In Fig. 1^ we plot the values of \eK\ versus the phase (^12 for 5i — 1 and the 
other parameters are chosen as in Fig. ^. From this figure, we conclude that the 
SUSY contribution with Hermitian Yukawa and universal soft scalar masses can not 
account for the experimentally observed indirect CP violation in the Kaon system. 
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In Ref. values for ek ^ 10"'^ have been obtained but these values require light 
gaugino mass (toi/2 ^ 100 GeV) which is now excluded by the new experimental 
limits on the mass of the lightets Higgs. Also it requires that the magnitude of 
the off-diagonal entries of the A-terms should be much larger (at least five times 
larger) than the diagonal ones, which looks unnatural. 

The chargino contribution to the — mixing is given by E3 

* \a.,b J i,j 

where Xi = (rn^+ /rriq)'^ . Here K is the CKM matrix, a, b are flavour indices, i, j label 
the chargino mass eigenstates, and V is one of the matrices used for djagonalizing 
the chargino mass matrix. The loop function H{xi,Xj) is given in RefE3. However, 
the chargino contribution can be significant only if there is a large LL mixing in 
the up- sector, namely Im((5]^^)2i ~ 10~^ and Rc((5|;^)2i ^ 10^^ E3. Such mixing 
can not be accommodated with the universal scalar masses assumption (i.e., 6i = 
1). In this case, the values of the Im((5]^^) are of order 10~^ which leads to a 
negligible chargino contribution to Ek- With non-universal soft scalar masses {Si ^ 
1) a possible enhancement in the chargino contribution is expected however, as we 
will show, this non-universality also leads to a larger enhancement in the gluino 
contribution. So, the dominant SUSY contribution to the K — K mixing in this 
class of models will be provided by the gluino exchange diagrams. 

A possible way to enhance Ek is to have non-universal soft squark masses at 
GUT scale E3. As mentioned above, the soft scalar masses are not necessarily uni- 
versal in generic SUSY models and their non-universality is not constrained by 
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Fig. 8. The value of \eii\ as a function of the parameters 5i (upper curve) and 82 (lower 
curve) for tan/3 = 5, and mg = = 200 GeV. 



the EDMs. Thus for 6i ^ 1 the mass insertion {5^2) RR is enhanced and we can 
easily saturate ex through the gluino contribution. To see this more explicitly, let 
us consider the LL and RR squark mass matrices in the super-CKM basis 

{^Drr-^'^ {Mlfy"- (63) 

Due to the universality assumption of Mq at GUT scale, the matrix ^AfJ^ 

remains approximately universal and the mass insertions (^12 )ll are sufficiently 
small (Im((5f2)LL ^ 10^^). However, since the masses of the squark singlets Af|, 
are not universal, Eq.(|4^), sizeable off-diagonal elements in are obtained. 

We find that Re(<5f2)_Rfi is enhanced from ~ 10"'^ in the universal case (5; — 1) 
to ^ for 5i ~ 0.7 while the imagenary part remains the same, of order 10^^. 
Thus, in this case, we have 

^\l^{{5i2)LL{5i2)RR) I ^ ^J\Re{5i2)RRlH5i2)LL\ ^ 10^^ 

which is the required value in order to saturate the observed result of Ek @- 

In Fig. ^ we show the dependence of \sk\ on the parameters 5i. There, the two 
curves, from top to bottom, correspond to the values of \£k\ versus 5i (for 82 = 1) 
and 82 (for 81 = 1) respectively. As explained above, in this scenario the main 
contribution of ek is due to LL and RR sectors and the LR sector has essentially 
no effect on £k ■ We also see that any non-universality between the soft scalar masses 
of the third and the first two generations can not lead to significant contribution 
to Ek and some spliting between the scalar masses of the first two generations is 
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necessary. This stems from the fact that the effect of the third generation on the 
mass insertion {5i2)rr is suppressed by V13 ~ 0(10^^) while V12 ~ sin^c- Finally, 
as we can see from this figure, in order to avoid over saturation of the experimental 
value of Ek, the parameter 61 should be of order 0.8. 



4.2. Direct CP violation 

Next let us consider SUSY contributions to e'/e. The ratio e'/e is a measure of 
direct CP violation in the K tttt decays and is given by 

e'/e = (imAo - - ImAa ) , (64) 

' V2|e|Re^o \ ^ J ^ ' 

where Aq,2 are the amplitudes for the A/ = 1/2, 3/2 transitions, and w = ReA2/ReAQ 
1/22 reflects the strong enhancement of AI = 1/2 transitions over those with 
AI = 3/2. Experimentally it has been found to be Re(e'/e) ~ 1.9 x 10""^ which 
provides firm evidence for the existence of direct CP violation. The ImAo.2 are cal- 
culated from the general low energy effective Hamiltonian for AS* = 1 transition, 

^eff^=' = EC4/^)0^ + /i.c. , (65) 



where Ci are the Wilson coefficients and the list of the relevant operators Oi for 
this transition is given in Ref.0@'S Let us recall here that these operators can be 
classified into three categories. The first category includes dimension six operators: 
Oi 2 which refer to the current-current operators. O-i-a for QCD penguin opera- 
tors and Oj-io for electroweak penguin operators E3. The second category includes 
dimension five operators: magnetic- and electricrdipole penguin operators Og and 
0.y which are induced by the gluino exchange Ea. The third category includes the 
only dimension four operator Oz generated by the sdZ vertex which is mediated 
by chargino exchanges E3. In addition, one should take into account Oi operators 
which are obtained from Oi by the exchange L R. 

In spite of the presence of this large number of operators that in principle can 
contribute to e'/e, it is remarkable that few of them can give significant contribu- 
tions. As we will discuss below, this is due to the suppression of the matrix elements 
and/or the associated Wilson coefficients of most of the operators. The SM con- 
tribution to e'/e is dominated by the operators Qe and Qg, and can be expressed 
as 

where — V*,Vid and the function F^r is given in Ref.El. By using our Hermitian 
Yukawa in Eq.(E^) we get 



e'/e ~ 7.5 x 10"''. 



(67) 
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Again, the SM prediction is below the observed value. Nevertheless, this estiniat 
can not be considered as a firm conclusion for a new physics beyond the SM since 
there are significant hadronic uncertainties are involved. 

The supersymmetric contribution to e'/e depends on the flavor structure of the 
SUSY model. It is known that, in a minimal flavor SUSY model, it is not possible 
to generate a sizeable contribution to e' je even if the SUSY phases are assumed 
to be large. Recently, it has been pointed out that with non-degenerate >l-terms 
the gluino contribution to e' j e can naturally be enhanced to saturate the observed 
value u. Indeed, in this scenario, the LR mass insertions can have large imaginary 
parts and the chromomagnetic operator Og gives the dominant contribution to e'/e 



^ , u , 11^3 ms F? 2 2 T 
Re — ^ „. ^, ,„ . mi mi Im 



^3 ~ 



(68) 



e J 647r2|£|ReAo rus + rud F^ 
where Cg is the Wilson coefficient associated with the operator Og, given by 

^9 = ^ i^i2)LL (^-^M^ix) - 3A/4(x)^ + ('5f2)Lfl^ (-^Mi(a;) " SAHx)^ 

where the functions Mi{x) can be found in Ref.i and X — mj./m'-. 

Using these relations, one finds that in order to saturate e'/e from the gluino 
contribution the imaginary parts of the LR mass insertions for a; ~ 1 should satisfy 
Im((5i2)Lfl ~ 10~^. Such values can easily be obtained in this class of models. How- 
ever, as mentioned above, in the case of Hermitian j4-terms and Yukawa couplings 
we have (<5f2)Lfl ^ iSi2)RL, hence we get Im Cg - Cg ~ Im [(^2)^^^ " (^12)^1^] - 
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10 ^ which leads to a negligible gluino contribution to e'/e E3. It is worth notic- 
ing that, due to the universality assumption of , the imaginary part of the 

mass insertion {Sf2)LL is of order 10~^. So its contribution to Cg is negligible with 
respect to the LR one which is enhanced by the ratio rrig/rag. To achieve the re- 
quired contribution to e'/e from the LL sector, one has to relax this universality 
assumptionnS to get Im((5f2)LL ~ 10^^. However, as we will discuss below, in this 
case the chargino contribution is also enhanced and becomes dominant. 

Now we turn to the chargino contributions. The dominant contribution is found 
to be due to the terms proportional to a single mass insertion e3. 

^""{jf =lui{^^K:2{Kb)LLK,}j F,,{xg^). (70) 

The function Fgi{xq^), where Xg^ — m-±/m?, is given in @. We find that the 
contributions involving a double mass insertion, like those arising from the super- 
symmetric effective sdZ, can not give any significant contribution, however we take 
them into account. The above contribution is dominated by (<5"2)l-L and in order 
to account for e'/e entirely from the chargino exchange the up sector has to employ 
a large LL mixing. Again, with universal M? lm{Si2)LL ^ 10"^ and the chargino 
contributions (as the gluino one) to e'/e is negligible. 
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Finally, we consider another possibility proposed by Kagan and Neubert to 
obtain a large contribution to e' j e c3. It is important to note that in the previous 
mechanisms to generate e' je one is tacitly assuming that A/ =1/2 transitions are 
dorninant and that the A/ = 3/2 ones are suppressed as in the SM. However, in 
RefJia it was shown that it is possible to generate a large e' je from the A/ = 3/2 
penguin operators. This mechanism relies on the LL mass insertion and 
requires isospin violation in the squark masses [jriu 7^ mj). In this case, the relevant 
AS" = 1 gluino box diagrams lead to the effective Hamiltonian 



J^cS 



V26 



D 

Y,{C^i^^)Q^+C^Q^) +h.C. (71) 
i=3 



where 



Ql — idaSa)v-A (9/39/3)y+A, Q2 — idaSfj)v~A iqf3Sa)v+A, (72) 
Q3 — idaSa)v-A i(lf3qp)v-A, Qi — idaSf})v-A (9/39a)y-A, (73) 

and the operators Qi are obtained fr'om Qi by exchaneing L ^ i?. It turns out that 
the SUSY AJ = 3/2 contribution to e'/e is given by 

Re ( 1 )^'^'^' ^ 19.2 f^^Vsf Im(<5f,),,. (74) 



Here, x^''^ = 
is given by 



e J \ nrig J 

2 rr. /rnr \2 



(^^) ' = (-Sj^) ' ^s'^ (™-) - 1 function K{x, y, z) 



32 2 

K{x, y,z)^ — [f{x, y) - f{x, z)] + — [g{x, y) - g{x, z)] , (75) 

where f{x,y) and g{x,y) are given in Ref.0. It is clear from the above equation 
that for mj^ = m^^, the function K{x^,x^,x^) vanishes identically. Thus, a 
mass splitting between the right-handed squark mass is necessary to get large 
contributions to e'/e through this mechanism. Furthermore, the lm{Sf2)LL has to 
be of order O{l0^^ — 10"^) to saturate the observed value of e'/e. It is clear that, 
with universal -/Vf? this contribution can not give any significant value for e'/e. 
Now we relax the universality assumption of A/? at GUT scale to enhance the 

mass insertion ((5i2)ll and saturate the experimental value of e'/e El. As men- 
tioned above, the non-universality between the first two generation of is very 
constrained by AMk and £k- Therefore we just assume that the mass of third 
generation is given by dsuiQ while the masses of the first two generations remain 
universal and equal to mg. This deviation from universality provides enhancement 
to both Ek and e'/e. We have chosen the parameters Si so that the total contribu- 
tions of ek from chargino and gluino are consistent with the experimental limits. 

In Fig. ^ we present the different gluino and chargino contributions to the e'/e 
and also the total contribution versus the imaginary part of the mass insertion 
((5^2 As explained above, there are two sources for the gluino contributions to 
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Fig. 9. The e'/e contributions versus the Imaginary part of the mass insertion {5i2)ll- 

e'/e: the usual A/ =1/2 chromomagnetic dipole operator and the new A/ = 3/2 
penguin operators. Additionally, there are two sources for the chargino contribu- 
tion to e'/e: the usual gluon and electroweak penguin diagrams with single mass 
insertion and the contribution from the SUSY effective sdZ vertex. As can be seen 
from this figure, the dominant contribution to e'/e is due to the chargino exchange 
with one mass insertion. It turns out that the chargino contribution with two mass 
insertions is negligible. As expected the gluino contribution via the chromomag- 
netic operator is also negligible due to the severe cancellation between the LR and 
RL contributions. The contribution from the A/ = 3/2 operators does not lead to 
significant results for e'/e. It even becomes negative (opposite to the other contri- 
butions) for Im((5f2)LL > 2.5 x lO'^. 

From this figure we conclude that a non-universality among the soft scalar 
masses is necessary to get large values of e'/e and Ek- 

5. CP violation in the B system 

5.1. Flavom — dependent SUSY phases and CP asymmetry in 
B JI^K, 

Recent results from the B-factories have confirmed, for the first time, the existence 
of CP violation in the i?-meason decays. Jn particular, the measurements of the 
CP asymmetry in the Bd i^Kg decay c3 have verified that CP is significantly 
violated in the i3-sector. In Bd — Bd mixing, the flavor eigenstates are Bd — (bd) 
and Bd — (bd). 

It is customary to denote the corresponding mass eigenstates by Bh = pBd+qBd 
and -Bl = pBd^qBd where indices H and L refer to heavy and light mass eigenstates 
respectively, and p= [1 + es)/ V^O'+l^sj), q= (1 - es)/ ^2{l + |es|) where Eb 
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is the corresponding CP violating parameter in the Bd — Bd system, analogous of 
e in the kaon system c3. Then the strength of Bd — Bd mixing is described by the 
mass difference 

AMb,=Mb^-Mb,. (76) 

whose present experimental value is AMb^ = 0.484 ± 0.010 (ps)~^ El. 

The CP asymmetry of the Bd and Bd meson decay to the CP eigenstate tpKs 
is given by 

a,Ksit) - . 3in(A^,,i). (77) 

i'Ksv ) r(BO(i)^ v^s) + r(BO(t)^ V-?fs) 

The most recent measurements of this asymmetry are given by 

av,Ks = 0.59 ±0.14 ±0.05 (BaBar) , 

a^Ks = 0.99 ± 0.14 ± 0.06 (Belle) . (78) 

where the second and third numbers correspond to statistically and systematic 
errors respectively, and so the present world average is given by a^pus — 0.79 ± 12. 
These results show that there is a large CP asymmetry in the B meson system. 
This implies that either the CP is not an approximate symmetry irunature and 
that the CKM mechanism is the dominant source of CP violation c3 or CP is 
an approximate symmetry with large flavour structure beyond the standard CKM 
matrix c3. Generally, AMb^ and a^Ks can be calculated via 

AAfB, = 2|(B°|i?,^/=2|50)|, (79) 

a^Ks = sin2/3eff , and /3eff = ^arg(B°|i?,^^=2|^°) . (80) 

where H^^^ is the effective Hamiltonian responsible of the AB — 2 transitions. 
In the framework of the standard model (SM) , a^Ks can be easily related to one of 
the inner angles of the unitarity triangles and parametrized by the Vck m elements 
as follows 

a^M^=sin2/3, /? = arg (^-|^) , (81) 

In supersymmetric theories the effective Hamiltonian for AB — 2 transitions, 
can be generated, in addition to the W box diagrams of SM, through other box dia- 
grams mediated by charged Higgs, neutralino, gluino, and chargino exchanges. The 
Higgs contributions are suppressed by the quark masses and can be neglected. The 
neutralino exchange diagrams are also very suppressed compared to the gluino and 
chargino ones, due to their electroweak neutral couplings to fermion and sfermions. 
Thus, the dominant SUSY contributions to the off diagonal entry in the B- meson 
mass matrix, 7Wi2(Bd) = {B°\H^^^=^\B'>^) , is given by 

Mi2{Bd) = Mlf{Bd)+MUBd)+MlUBd). (82) 
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where Mf^^Bd), Ml2iBd), and M^2 indicate the SM, gluino, and chargino cor 
tributions respectively. The SM contributioriis known at NLO accuracy in QCD t 
(as well as the leading SUSY contributions cHl) and it is given by 



Mf^iBd) - -^VBBBjl,MB,M^{VMfSo{xt), (83) 

where fs^ is the B meson decay constant, Bb^ is the renormalizatkm group invari- 
ant B parameter (for its definition and numerical value, see Ref. and reference 
therein) and rj — 0.55 ± 0.01. The function So{xt), connected to the AB = 2 box 
diagram with W exchange, is given by 

Axt - llx'^t+xf 3xf \nxt 
^'^'''>- A{l^x,r ~2il-x,r' ^^^^ 

where Xt = M^/M^. 

The effect of supersymmetry can be simply described by a dimensionless pa- 
rameter r? and a phase 29d defined as follows 



„2„2i 



Mi2iBd 



^,6-'" = (85) 



MffiBd)' 



where AMb^ = 2\Mi^{Bd)\r'^. Thus, in the presence of SUSY contributions, the 
CP asymmetry Bd — > V'-^s is modified, and now we have 

a^Ks = sin 2/3eff = sin(2/3 + 26*^) . (86) 

Therefore, the measurement of a^if^ would not determine sin 2/3 but rather sin 2/3off , 
where 



Mlf{Bd) 

: + 



2^. = arg(l+ -^/,,/„tM , (87) 



and MfysY(5^) ^ MUBd) + (Bd)- 

The most general effective Hamiltonian for AB = 2 processes, induced by gluino 
and chargino exchanges through AB ~ 2 box diagrams, can be expressed as 

5 3 

^eT=' = E c^{^^)Q^i^^) + E + ^-C , (88) 

i=l i=l 

where Ci(p), Ci{^) and Qi{^), Qiifi) are the Wilson coefficients and operators 
respectively renormalized at the scale /i, with 

Q4 = d^bt d{b% Qs = d%bl d^bl . (89) 

In addition, the operators Qi,2,z are obtained from (5i,2,3 by exchanging i ^ _R. In 
the case of the gluino exchange all the above operators give significant contributions 
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and the corresponding Wilson coefRcients are given by I 
CiiMs) = (24x/6(x) + mUx)) {6f,)i 



„2 



C2{Ms) - -,^^2MxUx){5%)hl, 

CsiMs) = -^36xUx)(df,)nL, (90) 
216m| 



RL 

C^iMs) = [[2Axfeix) + UOUix)) i6f,)LLiSf,)RR - mOMx)iSt,)LR{Sf,)RL 

where x — m^jm}- and is an average squark mass. The expression for the 
functions j^{x) and j>i{x) can be found in Ref.c3. The Wilson coefficients (71,2,3 are 
simply obtained by interchanging L <-> i? in the mass insertions appearing in Ci,2.3- 
Now we considerthe chargino contribution to the effective Hamiltonian in Eq. ( |88| ) 
in the mass insertion approximation. The dominant chargino exchange, can give sig- 
nificant contributions only to the operators Q\ and Qs in Eq.(^9|) lj. At the first 
order in the mass insertion approximation, the Wilson coefficients C'^^iMg) (by 
takin g d ifferent the mass of the stop-right from the average squark mass) take the 
form l3 

g?(Ms) ^ ^gg^',^, E {l^-il'l^^il' (('^^^)3i + 2A(^£j3i(^];j32) L-,{x.,,x,) 

- 2yt|y.i|V,ilA* ((^2L)3l('5]iL)31 + A(<5i;j32(<5]^j31 + A(J];j3l(5]^L)32)i?2(x„: 

+ Y,'VaV:,V,^V;, [{61^)1, +2X{SlMSRLh2 )r2{x^,x,,z)}, 

CUMs) = -^^^Y.Ur2Uj2V,^V., ((5EL)il+2A(<52j3l(<5EL)32)Lo(x,,X,), 

where Xi = m^+/?7i^, and the functions Lq{x, y) and L2{x, y) are given by 

Lo{x,y)^^^ ^y j' '"^^^ ^ (1 - .)^ (T^^ 

xh2{x) ~yh2{y) , 2 + 5x - Gxlna; 

i^2(x,,) = ' (l-.)3 +(13^ (93) 

and 

1 - , . 1 



R2{x,y,z) = (H^ix.z) - H2{y,z)) , R2{x,y,z) = (h2{x,z) - H2{y,z)) 

X — y X — y \ / 

H2{x,z) = — -U-l + x){x - z){-l + z){-l - X - z + 3xz) 
D2(x,z) L 
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-6 r 

H2{x,z) = <^ {-l + x){x - z)(-l + + (-2 + j;)z) 

D2{x,z) ^ 

+ 6a;2(-l + z)3log(x) -6(-l + a;)2z(-2a; + z + z2)log(z)| (94) 

Here D2(x, z) = (-1 + xf {x - z) (-1 + zf and i)2(a;, z) = (-1 + x)^ (a; - zf (-1 + zf . 
Notice that in the Hmit z I, both the functions R2{x, y, z) and R2{x, y, z) tend to 
L2{x,y). As in the gluino case, the corresponding results for Ci and Ci coefficients 
are simply obtained by interchanging L ^ R. in the mass insertions appearing in 
the expressions for Ci,3. 

In order to connect Ci{Ms) at SUSY scale Ms with the corresponding low en- 
ergy ones Ci{fi) (where /i ~ 0(mh)), one has to solve the renormalization group 
equations (RGE) for the Wilson coefficients corresponding to the effective Hamil- 
tonian in (p8|). Then, Ci(/i) will be related to Ci{Ms) by the following relationlla 

a(M) = EE (^^"^ + ^"^^■(^^)' (95) 

i s 

where Ms > rrit and 77 — as{Ms) / as{^j)- Thejalues of the coefficients 6-'^''*\ c[''^\ 
and ai appearing in ( p5| ) can be found in Ref.c2l. 

In tables (1) and (2) we present the bounds on real and imaginary parts of 
mass insertions respectively, by taking into account a light stop-right massc2l. We 
considered two representative cases of rhtj^ — 100, 200 GeV. Clearly, the light stop- 
right effect does not affect bounds on mass insertions containing LL interactions. 
From these results we could see that the effect of taking rhti^ < rh \s sizable. In 
particular, on the bounds of the mass insertions {i — 1,2) which 

are the most sensitive to a light stop-right. 

From the results in tables (10) and (11), it is remarkable to notice that, in the 
limit of very heavy squark masses but with fixed right stop and chargino masses, 
the bounds on the mass insertions ('^flx)3i('^fiL)3i ^^"^^ to constant values. This 
is indeed an interesting property which shows a particular non-decoupling effect 
of supcrsymmetry when two light right-stop run inside the diagrams in Fig. (1). 
This feature is related to the infrared singularity of the loop function R2{x, x, z) 
in the limit z ^ 0. In particular, we find that lim^^o ^^2(2;, a;, z) — f{x)/x, where 
X = m'^/rh'^, and f{x) is a non-singular and non-null function in x — 0. Then, in 
the limit rh >> the rescaling factor 1/m^ in will be canceled by the 1/x 
dependence in the loop function and replaced by l/m'^ times a constant factor. 

This is a quite interesting result, since it shows that in the case of light right 
stop and charginos masses, in comparison to the other squark masses, the SUSY 
contribution (mediated by charginos) to the AB = 2 processes might not decouple 
and could be sizable, provided that the mass insertions {6^1^)31 are large enough. 
This effect could be achieved, for instance, in supcrsymmetric models with non- 
universal soft breaking terms. 
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In the case of SUSY models with Hermitian flavor structure, we find that in most 
of the parameter space the chargino gives the dominant contribution to — 
mixing and the CP asymmetry ajj^Ks '^^^ the gluino is sub-dominant. As we 
emphasized above, in order to have a significant ghiino contribution for m ~ nig ~ 
500 GeV {i.e., toq ^ ^1/2 ^ 200 at GUT scale), the real and imaginary parts of 
mass insertion {5f.^)LL or {6f^)LR should be of order IQ-i and IQ-^ respectively. 
However, with the above hierarchical Yukawas we find that these mass insertions are 
two orders of magnitude below the required values so that the gluino contributions 
are very small. 

Concerning the chargino amplitude to the CP asymmetry aj/^n^, we find that 
the mass insertions (<53i)ijL and {S'^i)ll give the leading contribution to ajj^Ks- 
However, for the representative case of mo = mi/2 = 200 and ~ 7r/2 the values 
of these mass insertions are given by 

^|Imp-J3lP|=6xl0-^ (96) 

^|Im[(<5«J3lP|=4xl0-^ (97) 

^|Im[(5«^)3i(5»j32]| = 1 x 10-1 (98) 

These results show that, also for this class of models, SUSY contributions cannot 
give sizable eS'ects to aj/^i^g. As expected, with hierarchical Yukawa couplings 

(where the mixing between different generations is very small) , the SUSY contribu- 
tions to the B—B mixing and the CP asymmetry of i?d J/ipKs are sub-dominant 



31 



and the SM should give the dominant contribution. 

5.2. Flavom — dependent SUSY phases and CP asymmetry in 
B — > decays 

The Standard Model prediction for the CP asymmetry in rare B decays B Xs'^ 
is very small, less than 1%. Thus, the observation of sizeable asymmetry in this 
decay would be a clean signal of new physics. The most recentresult reported by 
CLEO collaboration for the CP asummetry in these decays is c3 

~ 9% < A'^p''' < 42% , (99) 

and it is expected that the measurements of A''^'^^ will be improved in the next 
few years at the i3-factories. 

Supersymmetric predictions for A''^''^ are strongly dependent on the flavour 
structure of the soft breaking terms. It was shown that in the universal case, as in 
the minimal supergravity models, the prediction of the asymmetry is less than 2%, 
since in this case the EDM of the electron and neutron constrain the SUSY CP- 
violating phases to be very small Furthermore, it is also known that in this case, 
one can not get any sizeable SUSY contribution to the CP-violating observables, e 
and e'/e. 

In this subsection, we explore the effect of these large flavour-dependent phases 
on inducing a direct CP violation in B Xg'^ decay. We will show that the values 
of the asymmetry A^*^ in this class of models are much larger than the SM 
prediction in a wide region of the parameter space allowed by experiments, namely 
the EDM experimental limits and the bounds on the branching ratio of _B ^ X^j. 
The enhancement of A^'''^ is due to the important contributions from gluino- 
mediated diagrams, in this scenario, in addition to the usual chargino and charged 
Higgs contributions. 

The relevant operators for & — 57 decay are given by 

g 

lOTT 



Qs = ■f^m,SLa^''T%nGl,. (100) 



El 

The expression for the asymmetry A^'*^ , corrected to next-to-leading order is 



given by @ 



2 



^_2z {viz) + b{z,6)) 



Im[C2C;] + Im[C7C;] 

z b{z,S)Im[C2Cl]y (101) 

where z = m^/m^. The functions v{z) and b(z,S) can be found in Ref.S. The 
parameter S is related to the experimental cut on the photon energy, Ey > (1 — 
S)mf,/2, which is assumed to be 0.9. 
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The SUSY contributions to the Wilson coefficients Cj^s are obtained by calcu- 
lating the b S7 and b ^ sg amplitudes at the electroweak scale respectively. The 
leading-order contributions to these amplitudes are given by the 1-loop magnetic- 
dipole and chromomagnetic dipole penguin diagrams respectively, mediated by 
charged Higgs boson, chargino, gluino, and neutralino exchanges. As pointed out 
above, SUSY models with non-universal A- terms may induce non-negligible con- 
tributions to the dipole operators Qy.s which have opposite chirality to Qt.s- In 
the MSSM these contributions are suppressed by terms of order 0{ms/mb) due to 
the universality of the A~terms. However, in our case we should take them into 
account. Denoting by Cr^g the Wilson coefficients multiplying the new operators 
Qt^s, the expression for the asymmetry in Eq.( |lOl| ) will be modified by making the 
replacement 

C,C* -> C.,C* + C.,C*. (102) 

The expressions for Cy^s are given in the appendix and C2 = (there is no 
operator similar to Q2 containing right-handed quark fields). Note that including 
these modifications ( 102 ) may enhance the branching ratio of i? — > Xg^ and reduce 
the CP asymmetry, since jCyp is replaced by jCyp -I- IC*?]^ in the denominator of 
Eq.(lOl). If so, neglecting this contribution could lead to an incorrect conclusion. 

In the EDM-free models we are considering, we found that the flavor dependent 
phase (/)23 gives a large contribution to the CP asymmetry. This casn simply be 
understood by using the mass insertion; the gluino contributions to C7 and Cg (and 
also Ci and C^) are proportional to ((^23) Li?, which receives a dominant contributions 
from A23 entry. The effect of the other flavor dependent phases on A^^p'^ is found to 
be very small. In Fig. |l^ we show the dependence of ^^^'^ on the phase (/)23 in the 
case of hermitian A-terms. We assume tan (3 — 10, the diagonal elements An = mo, 
and consider three values for the off-diagonal elements Aij,{i ^ j): \Aij\ = mo 
(curve 1), \Aij\ — 3mo (curve 2), and \Aij\ = 5mo (curve 3). The conditions of the 
correct branching ratio has been automatically imposed. 

We see that the CP asymmetry A^T^^'' can be as large as 15%, which can be 
accessible at the B facrories. Also as emphasised in Ref.^, larger tan/3 is, the larger 
the CP asymmetry A^^"'' become. In the case of symmetric A-terms, the CP phase 
in (^23)ifl is due to the CKM mixing only, however, it implies a considerable CP 
asymmetry, A^p"'' - 0(10%) 
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